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^^ ' Abstract: Bolometers have proven to be very good detectors to search for rare processes thanks 

■' ■ to their excellent energy resolution and their low intrinsic background. Further active background 

O ■ rejection can be obtained by the simultaneous readout of the heat and light signals produced by 

^^ particles interacting in scintillating bolometers, as proposed by the LUCIFER experiment. In this 

framework, the choice of the light detector and the optimization of its working conditions play a 

crucial role. In this paper, we report a study of the performances of a Germanium bolometric light 

detector in terms of signal amplitude, energy resolution and signal time development. The impact 



X 

_ ^_' of various operational parameters on the detector performances is discussed. 
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1. Introduction 

Bolometers [|I], |2|] have proven to be excellent detectors for application to the search for rare events, 
Uke neutrinoless double beta decay and dark matter interactions. This is mainly due to their good 
energy resolution and to the wide choice offered in terms of the absorber material. Using bolome- 
ters, it is often possible to match the characteristics of the detector material with the requirements 
of rare event searches. 

The strongest requirement, when coming to neutrinoless double beta decay (Ovj3j3), is the 
reduction of the radioactive background coming both from environmental sources and from con- 
taminations in the detectors and the surrounding materials [^]. In particular, a background assumes 
a crucial role [0|. The expected Ovj3j3 signal consists of two simultaneous electrons with a summed 
energy of a few MeV, an energy which can be released also by degraded a particles emitted from 
radioactive contaminations on the surface of the detectors or of the facing materials. The objective 
of the future experiments is challenging: the background level in the region of interest for Ovj3j3 
must be below 10^^ counts/(kg keV y). This cannot be achieved with standard bolometers, due to 
the lack of active a background rejection methods. 

Bolometers are sensitive calorimeters operated at ~ 10 mK where the detected signal is the 
temperature rise produced by the energy deposited in particle interactions. The various kind of 
particles (jS/ys, as, neutrons and recoiling nuclei) interact in different ways inside the absorber 
material but, in the hypothesis that any kind of energy deposition is eventually converted into heat, 
the heat signal recorded in the bolometer shows only a slight dependence on the particle type. 

Particle identification in bolometers can be obtained by the simultaneous readout of the heat 
and the light signal produced by scintillation in scintillating materials [^, ||, ^, |8[ ^, 10] or by 



Cerenkov effect in non-scintillating crystals [11]. Indeed, it has been observed in both cases that 



the amount of emitted light has a strong dependence on the particle type. 

The light detectors used so far for scintillating bolometers and for Cerenkov light readout 
consist in thin germanium slabs, also operated as bolometers. This choice was driven by several 
considerations. Being opaque semiconductors, they are sensitive over an extremely wide range of 
photon wavelengths and satisfy the very stringent radiopurity requirements of rare event searches. 
Their overall quantum efficiency can be as good as the one of photodiodes, providing at the same 
time an higher energy resolution. Moreover, they are much easier to operate at cryogenic tempera- 
tures. 

The operational characteristics of the light detectors (response function, energy resolution, 
noise, etc..) are crucial for the optimization of their overall performances, and consequently for 
the improvement of the background rejection capability. 

In this paper we present the results of the characterization of a light detector used for bolo- 
metric applications. The results were obtained in low temperature experimental runs carried out at 
Laboratori Nazionali del Gran Sasso (LNGS), near L'Aquila - Italy. 

2. Light detector operation 

We operated as bolometer a light detector (LD) produced by UMICORE, constituted by a disk- 
shaped pure Ge crystal (jZfSO mm x 300 /^m) grown using Czochralski technique. One side of the 
wafer was etched, the other one was polished. 

The temperature sensor of the LD is a3xl. 5x0.4 mm^ Neutron Transmutation Doped (NTD) 



Ge thermistor []12|], thermally coupled to the etched side of the crystal via 6 epoxidic glue spots of 
w 600 /im diameter and ^ 50 /xm height. 

The scintillation photons, as well as other particles, interact with the LD and deposit their 
energy within its volume. The energy deposition is eventually converted into phonons, causing a 
temperature rise of the system and, in particular, of the NTD thermistor. 

The thermistor works as a resistive device that converts temperature variations into resistance 
variations. Its resistance can be expressed as: 

R{T)=RoexpC^^ (2.1) 

where RqjTq and 7 depend on the doping level and are determined experimentally (typical values 
are: Rq = 1 .15 Q., To = 3.35 K and 7= 1/2). 

To measure the resistance variation, a bias voltage Vbias is applied across a pair of load re- 
sistors Rl/2 in series with the thermistor, as shown in Fig. |l]. The load resistance Rl is chosen 
much higher than the thermistor resistance R{T) for noise reasons. The resulting current in the 
circuit Ihoi remains almost constant when the temperature changes. The steady current produces 
a power dissipation which increases the temperature of the sensor and acts back on its resistance 



R{T), modifying its value according to eq. 2A. Any energy release in the wafer causes a further 
temperature rise. Since the thermistor is thermally coupled to the LD, this will result in a variation 
of R{T) and, therefore, of the voltage across the thermistor Vtoi, leading to a signal. Load resis- 
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Figure 1. Electric scheme of the bias circuit used for thermistor readout. 



tances used in this work range from 2 GD. to 1 1 GD., while the bias voltage was varied between 0.5 
and 30 V. 

The signals are amplified by means of Silicon JFETs and fed into an 18 bit NI-6284 PXI ADC 



unit [ |13| , |14| ]. The ADC is preceded by an analogic 6-pole low -pass Bessel-Thomson filter, with 
the roll-off rate of 120 dB /decade, to prevent aliasing effects on the acquired signal. The trigger 
is software generated on each bolometer and, when it fires, waveforms 250 ms long, sampled at 
8 kHz, are saved on disk. To maximize the signal to noise ratio, waveforms are processed offline 
with the optimum filter algorithm [ 
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2.1 Noise contributions 

In a scintillating bolometer, only a few percentage of the particle energy is converted into scintilla- 
tion photons. As a consequence, the corresponding light signals produced in the LD have energies 
of the order of few tens of keV [^, |l7]|. Dealing with small signals, it is fundamental to identify the 
different noise contributions affecting the measurements. 

The dominant noise contribution arises from mechanical vibrations of the cryogenic apparatus, 
propagated to LD and wires. Vibrations cause temperature fluctuations in the volume of the LD. 
These instabilities are a very dangerous source of noise since they have a frequency spectrum 
similar to the signal one, that is within the range — 100 Hz. In addition, the wire vibrations are a 
source of microphonic noise. The thermodynamic fluctuations of the light detector give a negligible 
contribution to the noise compared to the aforementioned sources. 

Other noise contributions come from electronics. The first one is the parallel Johnson noise 
of the load resistors, that develops across the dynamic impedance of the thermistor. As a conse- 
quence, its effect is larger at high values and becomes negligible at small values of the thermistor 
impedance. Considering the temperature of the resistors (room temperature) and the typical pa- 
rameters of our setup, the resulting contribution of this kind of noise is expected to be of the order 
of 1 — 30 eV rms. Further noise contributions are given by the preamplifiers and consist in a series 
white noise (30 — 120 eV rms) generated by the JFET resistances, and a series 1 // noise (5 — 25 eV 
rms). Being independent from thermistor impedance, these noises have negligible effects at large 
values and become dominant at small values of the thermistor impedance. The ADC quantization 
noise contribution was found to be minimal. We must remember that all the cited sources must be 
quadratically summed. 
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3. Experimental setup 

The characterization of the LD was structured in a set of measurements carried out in an Oxford 
200 ^He/^He dilution refrigerator located deep underground in the Laboratori NazionaU del Gran 
Sasso. The working temperature of the cryostat was varied between 10 and 20 mK. 

The LD was mounted in a copper frame and held in place by PTFE clamps, as shown in Fig. ^. 
The frame was accommodated in a container, also made of copper, connected to the refrigerator by 
means of a spring to reduce vibrations. 

To allow a proper calibration of the LD signal, a ^^Fe source, producing two X-rays at 5.9 and 
6.5 keV, was faced to the LD. The choice of such low energy calibration lines is due to the fact that, 
as previously mentioned, light signals produced in scintillating bolometers have typical energies of 
the order of 10 keV. 




Figure 2. Mounting of the LD into the copper frame. 



4. Detector performances 

We collected about 500 hours of data. The data-taking was split in short-lasting runs of few hours. 
As the main goal of the measurements was to characterize the LD and to optimize its performances, 
during these runs we systematically changed the working conditions, varying the bias voltage, 
the load resistance, the Bessel-Thomson cutoff frequency and the temperature. In the following 
sections we present the obtained results, trying to underline the effect of each operational parameter 
on the LD performances. 

4.1 Variation of the bias current 

Applying a Bessel-Thomson 3dB cutoff frequency fs = 200 Hz, we varied the bias voltage in the 
range 1 — 30 V, using two different load resistances: 2 and 1 1 G^. Acting on the polarization cir- 
cuit in this way, we were able to span a large range of bias current values, namely 150 — 15000 pA. 
We restricted our analysis to the events due to ^^Fe decay. An example of the Fe double peak, as it 
appears in the amplitude spectrum, is shown in Fig. ^ We performed a fit with a double Gaussian 
function, fixing the ratio of the peak positions and amplitudes to the known values from nuclear 



data [18] 



From the prominent ^^Fe peak we computed the signal yield as the ratio between the measured 
voltage variation across the thermistor and the energy of the line. The signal yield is expressed in 
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Figure 3. Fit of the ^^Fe X-ray double peak, obtained by biasing the LD with a current 4^/ ~ 4 nA. The 
relevant parameters of the fit are shown in the legend. Signal amplitude has been divided by the amplification 
given by the electronic chain, to obtain the voltage variation across the thermistor 
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Figure 4. Light detector performances resulting from the variation of the current in the bias circuit, using 
load resistances of 1 1 Gii (black dots) and 2 Gi2 (open circles). Both the signal yield (on the top) and the 
FWHM resolution (on the bottom) were computed on the 5.9 keV ^^Fe peak. 



/iV/keV. We also calculated the resolution of our detector as the FWHM of the ^^Fe peak in eV. 
Results are shown in Fig. Q 

As we can see from the picture, the maximum of the signal yield does not coincide with 
the best energy resolution. This is mainly due to the non-linear behavior of the thermistor. Indeed, 
increasing the signal yield also enlarges the noise in a not easily predictable way. As a consequence 
of the optimum filter application, the energy resolution depends on both noise power spectrum and 
signal shape, representing the true expression of the signal to noise ratio. Therefore, the optimal 
working point of the LD should be chosen as the one corresponding to the minimum of this curve. 
The value of the energy resolution of our detector at the optimal working point (that is not a single 
point but a large plateau between 2000 and 7000 pA) is ~ 250 eV. 

In both the plots of Fig. ^ we do not observe deviations from the common trends correspond- 
ing to variations of the load resistance, so we conclude that the value of the load resistance does 
not affect the signal amplitude and the energy resolution at the present level of accuracy. This is 
consistent with the a priori estimation of the load resistance Johnson noise discussed in sec. 0, since 
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Figure 5. Light detector noise power spectra for low (240 pA) and high (15000 pA) current, respectively in 
gray dotted line and in black solid line. Both the spectra were obtained averaging randomly triggered events 
including no pulses. 
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Figure 6. Light detector average pulse for low (240 pA) and high (15000 pA) current, respectively in gray 
dotted line and in black solid line. The reported waveforms were obtained averaging pulses belonging to the 
^^Fe lines. 



a contribution of ~ 1 — 30 eV rms is negligible with respect to the measured values of Fig. 0. The 
general decreasing behavior of the energy resolution is interrupted at very high polarization cur- 
rents. At these values, in fact, the parallel noise contributions become smaller than the ones given 
by the readout circuit (in particular, by the preamplifiers) and, from this point on, noise assumes the 
constant value of ~ 50 nV rms. Simultaneously, signal yield continues to decrease, determining the 
energy resolution worsening reported in the plot. It is important to stress that, however, a reduction 
of preamplifiers noise would not produce a sensible improvement of the resolution, since the signal 
yield decreases faster than the bias-related noise. 

Figures |5] and ^ show the noise power spectrum and the average detector response for low 
(240 pA) and high (15000 pA) bias current. The noise power spectrum was computed averaging 
randomly triggered waveforms including no pulses. The average pulse was calculated as the mean 
of a large number of pulses belonging to the ^^Fe lines. When comparing the noise power spectra 
we should keep in mind that the relevant part is the low frequency region (0 — 100 Hz) where the 
signal spectrum is confined. In this region, the increase of the bias current reflects into a reduction 
of the noise power spectrum of about two orders of magnitude. This is consistent with the energy 



Figure 7. Light detector signal rise time (black dots) as a function of the current in the bias circuit. The 
reported values were computed on pulses belonging to the ^^Fe lines. The trend is well reproduced by the 
convolution (black dashed line) of two time responses: one due to the 200 Hz Bessel-Thomson filter (gray 
dashed hne) and the other given by the readout circuit (gray dotted hne). 



resolution trend commented above. 

For what concerns the average detector response, the faster signal development in the case of 
higher bias is clear by looking at Fig. ^. We analyzed the distributions of the rise and decay times 
of ^^Fe pulses. 

The rise time is computed as the time difference between the 90% and the 10% of the leading 
edge. The mean value of the pulse rise time as a function of the bias current is shown in Fig. ^ 
The detector is remarkably faster (~ 3 times) at high bias currents. Speed response is limited by 
four different components. The first is introduced by the presence of the Bessel-Thomson filter, 
which gives a contribution that is independent from the thermistor impedance and, consequently, 
from the bias current. We evaluated this contribution by convolving in the frequency domain a step 
function with the transfer function of our Bessel-Thomson filter. This is reported in Fig. ^ by a 
dashed line. The second component is the rise time of the readout circuit, whose elements are the 
thermistor with its dynamic impedance, and the readout wires, having a parasitic capacitance of 
~ 250 pF and a negligible resistance. While the capacitance of the readout wires is constant, the 
working resistance of the thermistor depends on the bias current, as explained in section ^. The 
rise time of the readout circuit as a function of the bias current can be seen in Fig. ^ (gray dotted 
line). The last two contributions are given by the intrinsic thermal response of the absorber and the 
thermistor. They cannot be reliably estimated due to the high complexity of the phonon dynamics 
in the two components. Neveitheless the behavior of the experimental data is well reproduced 
by the convolution of the first two time responses (black dashed line in Fig. ^. Indeed, for low 
bias currents the RC contribution is dominant, while for higher currents the Bessel-Thomson filter 
accounts alone for the observed rise time. This implies that the thermal rise time, whatever it is, 
is negligible. When considering the rise time, it is important to underline that having fast signals 
reduces the probability of pile-ups that could spoil the detection efficiency. Therefore, large bias 
currents improve light detector performances in terms of both signal to noise ratio and detection 
efficiency. 

The decay time is computed as the time difference between the 30% and 90% of the pulse 
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Figure 8. Light detector signal decay time as a function of the current in the bias circuit (black dots). The 



evaluated decay times of the thermal pulses, evaluated by fitting the pulses with Eq. i.l as fit function, are 
also shown as open diamonds. All the reported values were computed on pulses belonging to the ^^Fe lines. 



trailing edge. The mean value of the pulse decay time as a function of the bias current is shown in 
Fig. |8[ As for the lise time, we can identify four different contributions to the decay time: Bessel- 
Thomson, readout, thermal response of the absorber and the thermistor. The first two have been 
already estimated for the rise time and are of the order of 2 — 6 ms, well below the decay time 
measured values. We deduce that the absorber thermal component is dominant in the decay time, 
provided that the thermistor contribution can be considered negligible. To confirm this hypothesis 
we performed a fit to the average detector response for a given bias current, in order to extract the 
decay time of the thermal component alone. The fit function is based on the thermal model response 
described in [19, ^. The fit accepts as input parameters measurable quantities like the sampling 
frequency, the gain, the Bessel-Thomson cutoff frequency, the bias voltage, the baseline voltage, 
the working resistance, the load resistance and the capacitance and gives as results the parameters 
of the thermal pulse which represents the temperature variation of the bolometer generated by an 
energy deposition: 

AT{t) = Ari-e^i^ + ae^^ {I - a)e'^) (4.1) 



Formula ^^ indicates that the thermistor temperature increases with one rise time, t^, and decreases 
with two decay constants, T^i and Tj2- The parameter a weights the two exponential decays and 
satisfies the condition < a < 1. The amplitude of the thermal pulse. At is also a function of the 
thermal elements and is directly proportional to the energy release E. An example of thermal pulse 
can be seen in Fig. ^ 

The decay time of the thermal component as a function of the bias current can be seen in Fig. || 
(open diamonds) and, as expected, it accounts entirely for the observed behavior. 



4.2 Variation of the Bessel-Thomson cutoff frequency 

Data were acquired using two different values of the Bessel-Thomson cutoff frequency: 200 Hz and 
550 Hz. Both the frequencies are well below the Nyquist limit of the sampling frequency (4 kHz), 
avoiding aliasing phenomena. The FWHM energy resolution as a function of the bias current is 
shown in Fig. fl^. The 200 Hz data are a subset of the ones already shown in Fig H. 
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Figure 9. Example of thermal pulse development for a bias current 4^/ = 150 pA. The amplitude is ex- 
pressed in the adimensional form -^P- (T is the LD working temperature). 
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Figure 10. Light detector FWHM resolution as a function of the current in the bias circuit. Data were 
acquired using two different Bessel-Thomson cutoff frequencies: 200 Hz (black dots) and 550 Hz (open 
circles). The reported values were computed on pulses belonging to the ^^Fe lines. 



We note that an higher Bessel-Thomson cutoff frequency causes a general worsening of the 
energy resolution, even if, applying a bias current around 3000 pA, we obtain the same minimum 
value of FWHM ~ 250 eV discussed above. 

For the same sets of data we have analyzed the distributions of the rise and decay times of 
^^Fe pulses. The mean values of the rise time as a function of the bias current is shown in Fig. 1 1 . 
As expected, a higher Bessel-Thomson cutoff frequency corresponds to a slightly faster rise time, 
especially at high bias currents. However, even considering the positive effect of a faster signal 
on the pile-up rejection, the general worsening of the energy resolution leads us to prefer lower 
Bessel-Thomson cutoff frequencies as the standard choice. 

No effect was observed on the decay time distribution where, as explained in the previous 
section, the thermal component is dominant. 



4.3 Variation of the working temperature 

In order to investigate the effect of the working temperature on the LD performances, we carried 
out several sets of measurements with various Ge light detectors, including the one presented so far 




Figure 11. Light detector signal rise time as a function of the current in the bias circuit. Data were ac- 
quired using two different Bessel-Thomson cutoff frequencies: 200 Hz (black dots, same data reported in 
Fig. 0) and 550 Hz (open circles). The reported values were computed on pulses belonging to the ^^Fe lines. 
Estimated contributions of 200 Hz(dashed hne) and 550 Hz (dotted line) Bessel-Thomson filters are also 
shown. 



*.*' 



800 

,[PA] 



Figure 12. Light detector FWHM resolution as a function of the current in the bias circuit. Data were 
acquired at two different cryostat temperatures: ^ 20 mK (black dots) and ^ 10 mK (open circles). The 
reported values were computed on pulses belonging to the ^^Fe lines. 



in this work, varying the working temperature of the cryostat from 20 mK to a significantly lower 
value. Among these measurements, we present here the one with the lowest achieved temperature 
(r^ 10 mK). The light detector used in this run is very similar to the one presented above, namely 
a Ge disk with dimensions jZfSO mm x 266 /im. It was mounted and operated in exactly the same 
way and, also, equipped with the same ^^Fe source. In Fig. 12, the FWHM energy resolution 



as a function of the bias current is shown for the two different values of the cryostat working 
temperature: ~ 20 mK and ~ 10 mK. 

With the exception of the low bias regime (< 200 pA), the higher temperature measurements 
clearly show better performances (roughly a factor 2 on the FWHM). We also note that the best 
energy resolution achieved using this detector at ~ 20 mK is comparable with the one of the other 



LD (see sec. 4.1) 
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5. Conclusions 

We investigated the perfonnances of Ge bolometric light detectors in terms of signal amplitude, 
energy resolution and signal time development. Data show a clear dependence on the applied 
bias current, favoring values around 2—7 nA. At higher polarization currents, series noise due 
to preamplifiers becomes dominant and limits the signal to noise ratio. The value of the load 
resistance has no influence on the signal amplitude and the energy resolution, confirming that the 
Johnson noise does not remarkably affect the light detector performances. Different choices of 
the Bessel-Thomson filter cutoff frequency do not critically change the LD behavior, provided 
that, however, too low or too high values would spoil the energy resolution. Temperature plays 
an important role: LD performances at ~ 20 mK are significantly better with respect to the ones 
observed at ~ 10 mK, the standard working temperature of bolometers for rare events searches. 

Acknowledgments 

The results reported here have been obtained in the framework of the Lucifer Experiment funded 
from the European Research Council under the European Unions Seventh Framework Programme 
(FP7/2007-2013) / ERC grant agreement n. 247115. Thanks are due to E. Tatananni, A. Rotilio, 
A. Corsi, B. Romualdi and F. De Amicis for continuous and constructive help in the overall setup 
construction. Finally, we are especially grateful to M. Perego and M. Guetti for their invaluable 
help. 

References 

[1] B. C. N. Booth and E. Fiorini, Low temperature particle detectors, Annu. Rev. Nucl. Part. Sci. 46 
(1996)471. 

[2] C. Enss and D. McCammon, Physical principles of low temperature detectors: Ultimate performance 
limits and current detector capabilities, J. Low Temp. Phys. 151 (2008) 5. 

[3] F. T. Avignone, S. R. Elliott, and J. Engel, Double Beta Decay, Majorana Neutrinos, and Neutrino 



Mass, Rev. Mod. Phys. 80 (2008) 481, [arXiv :0708.1033] 



[4] M. Clemenza, C. Maiano, L. Pattavina, and E. Previtali, Radon-induced surface contaminations in 
low background experiments, Eur Phys. J. C71 (2011) 1805. 

[5] P. De Marcillac et. al.. Experimental detection of a-particles from the radioactive decay of natural 
bismuth. Nature 422 (2003) 876. 

[6] C. Cozzini et. al.. Detection of the natural a decay of tungsten, Phys. Rev. C 70 (2004) 064606. 

[7] C. Arnaboldi et. al., CdWO^ scintillating bolometer for Double Beta Decay: Light and Heat 
anticorrelation, light yield and quenching factors, Astropart. Phys. 34 (2010) 143, 



[arXiv:1005.1239]. 



[8] J. Beeman et. al.. Performances of a large mass ZnSe bolometer to search for rare events, submitted to 



JINST [ arXiv:1303.4080 ] 



[9] J. Beeman et. al.. Performances of a large mass ZnMo04 scintillating bolometer for a next generation 
OvDBD experiment, Eur Phys. J. C 72 (2012) 2142. 



11- 



[10] J. Beeman et. al.. New experimental limits on the a decays of lead isotopes, Eur. Phys. J. A 49 (2013), 
no. 4 1-7. 

[11] J. Beeman et. al.. Discrimination of alpha and beta/gamma interactions in a Te02 bolometer. 



Astropart. Phys. 35 (2012) 558-562, [ arXiv : 1 1 6 . 62 8 6] 



[12] K. M. Itoh et. al.. Hopping conduction and metal-insulator transition in isotopically enriched 
neutron-transmutation-doped Ge:Ga, Phys. Rev. Lett. 77 (1996)4058. 

[13] C. Arnaboldi, G. Pessina, and S. Pirro, The cold preamplifier set-up ofCUORICINO: Towards 1000 
channels, Nucl. Instrum. Meth. A 559 (2006) 826. 

[14] C. Arnaboldi et. al.. The front-end readout for CUORICINO, an array of macro-bolometers and 
MIBETA, an array ofmu-bolometers, Nucl. Instrum. Meth. A 520 (2004) 578. 

[15] E. Gatti and P. F. Manfredi, Processing the signals from solid state detectors in elementary particle 
physics, Riv. Nuovo Cimento 9 (1986) 1. 

[16] V. Radeka and N. Karlovac, Least-square-error amplitude measurement of pulse signals in presence 
of noise, Nucl. Instrum. Methods 52 (1967) 86. 

[17] C. Arnaboldi et. al.. Characterization ofZnSe scintillating bolometers for Double Beta Decay, 



Astropart. Phys. 34 (2011) 344, [arXiv: 1006 . 2721 1. 



[18] S. Y. F. Chu, L. P. Ekstrom, and R. B. Firestone, The Lund/LBNL nuclear data search. Available at 
nucleardata. nuclear lu. se/nucleardata/toi. 

[19] M. Carrettoni and M. Vignati, Signal and noise simulation ofCUORE bolometric detectors, JINST6 



(201 1) P08007, [ arXiv :1106.3902] 



[20] M. Vignati, Model of the Response Function of Large Mass Bolometric Detectors, J. Appl. Phys. 108 



(2010)084903, [arXiv: 1006 . 4043). 



-12- 



